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ABSTRACT 

The  Haltiner  theory  of  saturated  steady  state 
convective  currents  which  Includes  the  transfer  of  heat 
and  momentum  by  lateral  diffusion  as  well  as  the 
entrainment  of  environmental  air  is  discussed.   A 
modified  model  is  presented  which  includes,  as  a 
separate  term,  the  drag  of  the  liquid  water  carried 
along  by  the  ascending  current .   The  system  of  equations 
is  integrated  numerically  for  several  cases  previously 
considered  by  Haltiner  allowing  entrainment  to  assume 
negative  values  (detrainment )  and  the  results  are 
compared „   The  results  are  in  good  agreement,  and  show 
that  the  net  effect  of  the  opposing  actions  of  liquid 
water  drag  and  detrainment  is  a  reduction  in  the 
computed  vertical  velocities.   This  reduction  appears 
to  be  a  correction  in  the  right  direction  on  the 
basis  of  a  comparison  to  some  observed  data. 

The  writer  wishes  to  express  his  appreciation  for 
the  assistance  and  encouragement  given  him  by 
Professor  GoJ.  Haltiner  of  the  U.S.  Naval  Postgraduate 
school  in  this  investigation. 
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lo    Introduction 

The  concept  of  simple  "adiabatic  parcel  ascent"  as 
a  first  approximation  to  convective  phenomena  has  long 
been  in  use;   it  has,  however,  a  number  of  shortcomings. 
As  a  result  of  evidence  (Stommel, [s]  )  that  a  rising 
convective  current  acts  like  a  jet,  dragging  in  surround- 
ing air  as  it  progresses  upward,  Austin  and  Fleisher 
suggested  that  this  process,  called  "dynamic  entrain- 
ment",  is  necessary  to  satisfy  mass  continuity.  C'] 
Haltiner [h]  summarized  the  work  of  Houghton  and  Cramer, C5"] 
Bunker, [2!  and  Malkus,f£T]  and  suggested  a  theory  of 
saturated  convective  currents  which  includes  the  effects 
of  lateral  diffusion  of  heat  and  momentum,  and 
extended  computations  to  elevations  typical  of  large 
cumulus  clouds c 

Haltiner  [Hj  derived  a  system  of  equations  which 
describe  a  saturated  convective  current.   The  derivation 
of  these  equations  is  essentially  the  same  as  the 
derivation  of  the  modified  system  of  equations  presented 
belowo   The  Haltiner  model  assumes  that  during 
detralnment  the  cross-sectional  area  of  the  rising 
column  remains  constant.   However,  when  mass  continuity 
considerations  calls  for  the  ejection  of  mass 
(detrainment )  s.  a  second  system  of  equations  is  applied 
based  on  the  assumptions  that  the  rate  of  entrainment 
remains  identically  zero  and  mass  continuity  is 


satisfied  by  an  increasing  cross-sectional  area  in 
the  rising  air„   In  this  paper  the  detrainment  phase 
will  be  treated  by  retaining  the  assumption  that  the 
cross-sectional  area  of  the  rising  air  remains  constant, 
thus  permitting  the  entrainment  to  assume  negative 
values o 

The  effect  of  the  drag  of  the  liquid  water  carried 
along  in  the  rising  air  was  not  explicitly  considered 
in  the  Haltiner  model,  although  it  can  be  qualitatively 
included  by  adjusting  the  rate  of  diffusion  of 
momentumo   Here  the  drag  effect  of  the  liquid  water  will 
be  specifically  included  as  a  separate  term. 

The  following  notation  will  be  used  throughout 
this  paper: 
T     temperature 
z     height 
q     specific  humidity 

cp    specific  heat  of  air  at  constant  pressure 
L     latent  heat  of  vaporization 
g     gravity 
w     velocity 
k     diffusion  coefficient 
Rd    gas  constant  for  dry  air 
R     gas  constant  for  water  vapor 
N     In  M 
M     Mass  rate  of  flow  per  unit  time 


I        "specific  humidity"  of  liquid  water 
J"   cross-sectional  area  of  rising  column 
/*       density 
p   pressure 
%a      lapse  rate 
subscripts: 
e   environment 
v   virtual 
o    initial  value 
2,    The  system  of  equations 

a.   The  continuity  equation. 
Consider  a  vertical  current  of  air  in  which  the  velocity, 
w,  is  horizontal!"  uniform.   If  the  field  of  vertical 
velocity  is  not  uniform  w  may  be  considered  to  be  the 
mean  value  in  the  horizontal  at  any  given  level.  The 
mass  of  air  per  unit  time  flowing  past  any  given  level,  z  , 
is  given  by  ^A 

where  M  is  the  mass,/'  the  density  and  0"  the  cross- 
sectional  area  of  the  vertical  current.   Under  steady 
state  conditions  the  mass  at  an  adjacent  level  Z  +dzis 
the  sum  of  the  mass  rate  of  flow  at  level  2  ,  and  air 
entrainment  fromi  to  z  +  dz  .   Thus, 


d  H     d/<>         cJ  w    ci  Q- 


Defining  a  parameter  N—  InM,  we  may  write 

(2)   ^=J_i^+-L^4--Lck 
d  z.  ,  °     J  z  w  dz.     (T  dz 


Differentiation  of  the  equation  of  state  p  -  Rj/*3""^ 
yields 

/'dz    p  dz    T  dz 

We  next  make  the  assumption  that  at  any  level  the 

pressure  in  the  rising  air  is  the  same  as  that  of  the 

environment,  then   _L  H£  -  _L  2)Q 

P  dZ    jo  cHi 

From  the  hydrostatic  equation  dp/£>z-:=  ~  °/%        it 

follows  that 


•  o\   p        |  ^j3  D 

P  dz      ft  3*        Rd'Xe 


We  may  now  eliminate  the  density  term  from  (2)  by 
substitution  from  (3)  and  (4)  to  yield  the  continuity 
equation  for  the  rising  air 

(5)   dN=_LcUr+J_ck-_  J _L  dX 

dz    w  dz  (7  dz  ^d~^    X  ^z 

bo   The  Equation  of  Motion. 

The  changing  mass  of  the  rising  air  due  to  entrainment 

or  detrainment  requires  the  use  of  an  equation  of 

motion  for  a  variable  mass.   If  a  mass  M  is  moving  at 

velocity  w  at  time  t  combines  with  a  second  mass  dM 

moving  at  velocity  w  during  a  time  interval  dt 


and  w  -»-  dw  represents  the  velocity  of  the  combined  masses 
at  time  t  -*-  dt,  then  the  momentums  at  time  t  and  t  +  dt 
are  related  by  the  equation 

Here  E  represents  the  external  forces  per  unit  mass. 
This  relationship  yields  the  equation  of  motion 

The  external  forces  acting  on  the  system  include  the 
pressure  and  gravity  forces,  which  may  be  expressed  as 

9(~C"""Xe)/  X  '  and  frictional  forces  F, 

the  latter  including  the  drag  of  the  liquid  water  carried 
by  the  ascending  column  of  air.   The  equation  of  motion 
for  the  liquid  water  rising  with  the  air  may  be  written 

(7)  idw  =  _  Jcj  +D 
d  t 

where  [    is  the  mass  of  liquid  water  per  unit  mass  of 
air  which  we  will  term  the  "specific  humidity"  of  the 
water  and  D  represents  the  drag  of  the  air  acting  on 
the  liquid  water  of  massj?  .   The  negative  of  D  is  the 
drag  of  the  liquid  water  on  a  unit  mass  of  air.   From 
eq  ( 7) ,  we  have 

substituting  in  eq  (6)  we  obtain 

)-  Jw 
M 


simplifying  we  have 


assuming  steady  state  conditions  we  obtain 

W  d  w        ' 


(8)     el*     (M) 


c,   Lateral  diffusion  of  heat  and  momentum 
In  addition  to  the  exchange  due  to  entrainment, 
turbulent  exchange  of  heat  and  momentum  must  also  be 
consideredo   Priestley  C7J ( 1956)  investigated  the 
mechanisms  of  diffusive  heat  and  momentum  exchange  in 
a  plume  of  hot  gases „   Haltiner  [*f'J  suggested  that, 
except  for  the  scale  of  the  phenomena,  the  diffusive 
exchange  processes  in  a  convective  current  are  similar 
to  the  processes  defined  by  Priestley.   The  effects  of 
lateral  diffusion  were  included  in  the  Haltiner  model 
by  the  extention  of  the  Priestley  equations  to  convective 
currents 0   The  heat  loss  per  unit  time  is  taken  to  be 

(9)    -M,cfx  (e-ej 

or  approximately 

(lO)-yfe,  Cp  (T-"t) 

It  is  also  necessary  to  saturate  the  environmental  air 
entrained  into  the  updraft  and  a  term  of  the  form 

(id  -Jk.u  (a-  4e) 

is  required^,  where  L  is  the  latent  heat  of  vaporization. 
The  loss  of  momentum  per  unit  time  per  unit  mass  is 
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similar  to  that  of  the  loss  of  heat  and  is  assumed  to 
be  proportional  to  the  difference  between  the  velocity 
of  the  convective  current  and  the  environment,  with  the 
latter  taken  to  be  zero,  with  the  result  —  J^w  where  H{ 
has  the  dimension  second".   With  this  form  for  the 
friction  term,  eq  (8)  becomes 

<«>  6fm  -  rfjtt[&(T^-*)-l*'-w«^  -K~ 

The  dimensions  of  k(  imply  that  diffusive  losses  are 
independent  of  updraft  velocity.   Haltiner  suggests  a 
second  approach  to  the  problem  of  diffusion  which  makes 
the  loss  of  heat  and  momentum  dependent  upon  updraft 
velocity0   The  first  approach,  however,  will  be  employed 
in  this  paper „ 

do   Thermodynamic  equation. 

In  describing  the  thermodynamic  process  undergone 
by  the  rising  air  a  number  or   factors  must  be  included. 
The  following  must  be  taken  into  account  in  accordance 
with  the  First  Law  of  Thermodynamics:  the  heat  released 
by  condensation,  the  heat  required  to  saturate  the 
entrained  environmental  air  and  to  raise  its  temperature 
to  that  of  the  rising  column,  the  heat  lost  to  the 
environment  through  diffusion,  the  internal  energy 
changes  and  the  work  done  during  expansion.   The 
resulting  equation  embodying  these  factors  is  as 
follows? 


(13)  -ML*H   -Uq-<fc)dM  -cr(T-Te)dM 

The  terms  on  the  left,  in  the  order  of  their 
occurrence,  represent  respectively,  (a)   the  heat 
released  in  condensation;  (b)   the  heat  required  to 
saturate  the  entrained  air,  (c)  the  heat  required  to 
bring  the  entrained  air  to  the  temperature  of  the  rising 
column,  (d)  the  heat  lost  through  diffusion  in  the 
interval  dt  required  for  the  air  to  rise  a  distance 
dz<,   The  right  side  of  the  equation,  represents  the 
combined  internal  energy  change  and  work  by  the  gas  in 
expanding  from  pressure  p  to  pressure  p-vdp.   By 
substitution  from  eq  (*0,  eq  (13)  becomes,  with 
appropriate  manipulation, 

-  Ld^  _  [c,(T-T.VUci-*>]d*! 

lw  az  dz 


J  2     Xe 

ec   Moisture  and  liquid  water. 

It  has  been  assumed  that  the  rising  air  is  sat- 
urated; the  specific  humidity  of  the  rising  air  is  then 
a  function  of  temperature  and  pressure.   As  a  good 
approximation  we  have  the  expression  <3  =  .bXl^/p     where  e 
is  the  saturation  vapor  pressure.   By  differentiation 
this  becomes 
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(15)  I  dq  __   i_  de  _  J.  le 
q  dz     c  cU  f  ciz 

The  last  term  in  eq  (15)  may  be  replaced  with  9^'^^ve 

from  eq  (9).  and  from  Clapeyron's  equation,  _  -  Ld^1. ,  . 

v 

where  R  is  the  gas  constant  for  water  vapor,  an 


alternate  expression  for  the  second  term  is  obtained, 
yielding 

U6)  qd*  "   R,T*J*     R,Tve 

The  liquid  water  content  of  the  cloud  is  the  amount  of 
moisture  condensed  less  the  amount  required  to  saturate 
the  environmental  air  entrained  and  that  exchanged 
through  diffusion.   Defining  X    as  the  "specific  humidity" 
of  the  liquid  water  we  may  write 

(17)  d(Ml)=-  rtd<i  -(<±-ctfdW  -MVj^^-cjpidt 

expanding  eq    (1?)    yields 


dz     dz  dz 


w 


Equation  (1?),  however,  is  not  valid  during  detrainment 
and  the  change  in  liquid  water  content  must  be  piece- 
wise  defined.   For  the  detrainment  phase  we  have 

(19)  d(M«  =  -M4*  +3!divl-  MVi,(<i-<k)dt 

whicV  upon  expansion  leads  to 


(20)   Sj       _    dq      K     (<l~(\e) 

Two  additional  relationships  between  temperature  and 
virtual  temperature  will  be  required.   One  of  these  is 
the  well  known  relationship 

(21)  T  =(l  +  .6l<0T 


the  second  relationship  is  obtained  from  eq  ( "1 )  by 
differentiating  with  respect  to  z ;  this  yields 


Substituting  from  eq  (16 )  this  may  be  written 

By  using  some  mean  values  this  may  be  approximated  by 


(22)   dT, 


/l  +  1  2.S~q\  cjT 


f,   The  environment. 

The  environmental  parameters  1^  ,  Tve   ,  qe ,  and  w 
appear  in  one  or  more  of  the  equations  derived  above; 
it  is  therefore  necessary  to  define  the  environment  in 
some  way0   We  have  wide  latitude  in  choosing  the  struc- 
ture of  the  environment;  it  may,  for  instance,  conform 
to  an  actual  sounding.   A  simple  model  was  used  in  this 
investigation.   The  lapse  rate  Xe  and  the  relative 
humidity  were  assumed  constant,  and  \       was  assumed  to 

10 


be  zero.   Under  these  assumptions  T^   may  be  expressed 
as 

(23)   Te  -  Teo-ye* 

The  environmental  specific  humidity  may  be  determined 
as  follows:  From  eq  (16)  the  saturation  specific 
humidity  is  given  by 

and  since  the  relative  humidity  is  constant  c^  =  p^,  we 
may  write 


(24)         ^   R  ^   ■ 

If  the  lapse  rate  &e  ,  the  initial  temperature  T^  ,  and 
the  relative  humidity  are  chosen,  then  the  values  of TG , 
"T   ,  and  qe  may  be  obtained  at  any  level  from  equations 
(21),  (23),  and  (24). 

3<>    The  complete  model. 

Having  defined  the  environment  in  sub  sbfct&dtoJD 
equations  (5),  (12),  (14),  (16),  (21)  and  (18)  or  (20) 
constitute  a  system  of  six  equations  in  7  unknowns,  w, 
T  ,'T  iq  .  i  >   N,  and  <J    .   We  must  now  establish  a 
relationship  defining  <T  .   The  simplest  solution,  and 
the  one  which  will  be  employed  here,  is  to  assume  that 

(T  is  identically  zero.   This  assumption  is  quite 
likely  an  oversimplification,  but  on  the  other  hand 
presently  available  information  on  the  behavior  of  the 
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cross-sectional  area  of  an  updraft  is  scanty  and  somewhat 
conf lictingo 

The  system  of  equations  (5),  (12),  (1*0,  (16),  (21), 
and  (18)  or  (20)  may  be  solved  first  for  the  cloud 
temperature  lapse 

ciT 


(25) 


dz 


and  then  successively  for  the  other  parameters 


(26)  d- 


4  ve 


(27) 


dw- 


CI+JP) 


»  -+- 


T.-T 


'vg 


*e_(' 


(28)  -p 


(29a) 


41 


RJa  dz    R^ 


_  ^  U-5-) 


w 


or 


(29b)  ^?=_  li  _M  q-<j.) 

cl?     d  2- 


w 


Equations  (21)  through  (24)  are  used  as  needed.   Given 
initial  values  of  T,  N,  w ,  $    ,  ~Q  ,  q  ,  c^ ,  the  foregoing 
system  of  equations  may  be  integrated  by  numerical 
methods  over  the  desired  range  of  z. 
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Since  adequate  data  concerning  the  vertical  velocity 
of  the  entrained  environmental  air  is  lacking;  eq  (2?) 
will  be  further  simplified  by  assuming  that  we  = o . 
Moreover  since  typical  values  of  J?  are  on  the  order  of 
o005  we  may  place  \/(\  +  i \)   =  1  .   The  error  introduced 
by  the  latter  approximation  will  be  on  the  order  of  \ 
of  one  percent o 
If  eq  (27)  is  modified  as  indicated  above, it  becomes 

(30)  ^T  =  _wdw  _  J?,  +  J.  /X-X«   _  i) 


c>  z 


a  form  which  differs  from  the  corresponding  equation 
in  the  Haltiner  model  only  in  the  last  term,  which  has 
the  form 
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1  •£ 


Similarly  in  the  Haltiner  equations  for  d~T/dz    and 
dlVdz,  the  term  2.  /v  ~~%A  is  replaced  by  -i-  /VTvej\ 

in  the  modified  model.   Equations  (28)  and  (29a)  are 
exactly  the  same  as  their  counterparts  in  the  Haltiner 
model o   The  only  other  difference  in  the  two  models  is 
the  appearance  of  the  additional  eq  (29b). 

^0    Computations  and  Results. 

The  foregoing  system  of  equations  was  integrated 
by  numerical  methods  on  an  electronic  digital  computer, 
the  National  Cash  Register  102A.   With  a  given  set  of 
initial  values,  properly  scaled  inputs  to  the  integration 
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routine  were  computed  for  the  model  given  above.   These 
values  were  then  integrated  over  an  interval  of  10  meters 
yielding  a  new  set  of  "initial"  values  for  the  level  z  + 
10  meters  which  served  as  an  input  for  the  computation 
of  a  new  set  of  inputs  for  the  integration  routine. 
Computation  proceeded  in  this  stepwise  fashion  until  the 
top  of  the  cloud  was  reached. 

Computations  were  made  for  four  sets  of  initial 
conditions.   In  each  case  the  initial  values  corresponded 
to  one  of  the  sets  of  initial  values  used  by  Haltiner^H^ 
In  Case  (a)  and  Case  (b)  the  same  set  of  initial  values 
were  used,  but  in  Case  (a)  only  the  effect  of  detrainment 
was  considered.   The  effect  of  the  drag  of  the  liquid 
water  was  eliminated  from  Case  (a)  by  an  appropriate 
modification  to  the  complete  program.   In  Cases  (b), 
(c),  and  (d),  both  the  effects  of  detrainment  and  the 
drag  of  the  liquid  water  are  Included. 

The  initial  values  used  in  the  computation  are  as 
follows: 
Case  (a) 

T0  s  20C,  Teo  "  19°C,  w0  = 1  meter/second 

k-j_  =   ,001/second,  ^  e  =.7°C/100  meters 
r  =   80$ 
Case  (b) 

Same  as  Case  (a)  except  that  both  detrainment  and 
water  drag  are  included. 
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Case  (c) 

Same  as  Case  (b)  except  that  \^     .6°  C/10C  meters 
Case  (d) 

Same  as  Case  (b)  except  that  r=  90$ 

The  results  of  the  computations  are  presented  in 
figures  1  through  6  and  table  1. 

Figures  1  and  2   present  the  results  of  Cases  (a) 
and  (b),  in  which  the  parameters  W,  #    ,  M/M0  and  A  T 
are  plotted  as  functions  of  height.   The  parameter  M/M0 
is  the  cumulative  change  in  mass  experienced  by  a  unit 
of  air  due  to  entrainmento   A  value  of  two,  for  example, 
indicates  that  the  mass  has  doubled.   The  parameter  A  T 
is  defined  as  T   -  Te  . 

Figures  3»  ^»  5»  and  6  present  a  comparison  of  the 
results  obtained  by  Haltiner  \^4]  f  or  the  set  of  values 
listed  as  Case  (a)  with  the  results  of  Cases  (a)  and  (b) 
Each  parameter  Is  plotted  separately  as  a  function  of 
height,.   The  values  computed  by  Haltiner  and  the  values 
from  Case  (a)  should  be  the  same  at  and  below  the  level 
where  detrainment  begins.   The  systematic  difference  in 
two  values  below  this  point  is  the  result  of  minor 
differences  in  the  computational  procedures  (primarily 
in  the  computation  of  the  environmental  parameters) 
followed  in  this  investigation  and  those  employed  by 
Haltiner,,   This  systematic  difference  is  not  significant 


15 


in  Itself,  but  it  does  illustrate  how  small  differences 
at  one  level  propagate  upward  and  tend  to  accumulate* 

The  effects  of  detrainment  and  the  drag  of  the 
liquid  water  are  best  illustrated  in  figures  k   and  5» 
It  is  apparent  that  the  effect  of  detrainment  alone  is 
to  accelerate  the  vertical  current.  The  cloud  height 
is  greater  in  Case  (a)  than  it  is  in  either  Case  (b)  or 
the  corresponding  case  computed  by  Haltiner\_H]  .   It  is 
also  apparent  from  the  figures  that  in  Case  (a)  the 
magnitude  of  the  maximum  values  of  w  and  M/M0  and  the 
elevation  at  which  these  maxima  occur  are  greater  than 
either  the  values  obtained  by  Haltiner  or  the  values 
of  Case  (b)0   These  results,  and  the  results  of  Cases 
(c)  and  (d)  are  summarized  in  table  one. 

5o    Conclusions 

The  number  of  cases  investigated  was  not  exhaustive, 
but  was  sufficient  to  show  that  the  modified  model 
yields  results  that  conform  quite  closely  to  the  results 
of  Haltiner[_4-]  •   In  particular,  the  cloud  heights  were 
substantially  the  same  but  the  vertical  velocity  and  the 
amount  of  entrainment  are  reduced.   These  results  are 
the  natural  consequence  of  the  combination  of  detrainment 
and  the  drag  of  the  liquid  water.   Detrainment  on  the 
one  hand  provides  an  acceleration,  while  the  drag  of  the 
liquid  water  acts  as  a  retarding  force.   These  forces 
tend  to  balance  one  another,  but  the  drag  of  the  liquid 
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water  somewhat  overba lances  the  acceleration  due  to 
de trainmen t „ 

Haltiner  V.  ^  J  investigated  an  actual  cloud  for  which 
observational  data  was  available,  and  found  that  his 
results  were  in  good  agreement  with  the  observed  values. 
He  did  find,  however,  that  his  values  of  vertical 
velocity,  although  of  the  right  order  of  magnitude,  were 
somewhat  in  excess  of  the  observed  values.   To  the 
extent  that  the  modified  model  reduces  vertical  velocity, 
the  modifications  appear  to  be  a  correction  in  the  right 
direction. 
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liquid  water  content  JJ  and  teniueratnre  exeoss  over  the 
environ  lent  A  T  as  a  function  oi  height  tor  Case  (a). 
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Figure    2.      Vertical    velocity   v;,    tro.il  t  i  onal    raass    increase  I'/:'e, 
liquid   water   content  Jj    and    tenperattire    excess    aver   the 
environ  rent  A  T   as    a  function  of  height  for  case  (b). 


19 


9< 

81 


1  *  ■ .  A  T  (°  C ) 


Figure  "J> .      Tempera  bur  f«  excess  over  eiivirou..ic.i  t  A  T  as  a 
function  of  height  for  Case  (a),  dotted,  Case  (u),  solid, 
and  after  llaltiuer,  dashed. 
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Figure  5«   Vertical  velocity. w,  as  a  function  of  height  for 
Case  (a),  Case  (b)  and  after  Haltiner.   Tick  marks  on 
curves  indicate  level  at  which  dctrainment  started. 
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Figure  6.   Liquid  water  content  as  a  function  of  height  for 
Case  (a),  dotted,  Cose  (b) ,  solid,  and  after  Haltiner,  dashed. 
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